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Abstract 


This research aims to optimize the mixing process in gas-lift anaerobic digesters of municipal sewage sludge 
since mixing and maintaining uniform contact between methanogenic bacteria and nutrients is essential. 
Wastewater municipal sludge sampling was performed at the Ahvaz West treatment plant (Chonibeh, Iran) 
during the summer of 2022. A Computational Fluid Dynamics (CFD) model was implemented to simulate, 
optimize, and confirm the simulation process using ANSYS Fluent software 19.0. The velocity of the inlet-gas 
into the digester was determined and a draft tube and a conical hanging baffle were added to the digester design. 
Different inlet-gas velocities were investigated to optimize the mixing in the digester. Furthermore, turbulence 
kinetic energy and other evaluation indexes related to the sludge particles such as their velocity, velocity 
gradient, and eddy viscosity were studied. The optimal inlet-gas velocity was determined to be 0.3 ms”. The 
simulation results were validated using the Particle Image Velocimetry (PIV) method and the correlation 
between CFD and PIV contours was statistically sufficient (98.8% at the bottom corner of the digester’s wall). 
The results showed that the model used for simulating, optimizing, and verifying the simulation process is valid. 
It can be recommended for gas-lift anaerobic digesters with the following specifications: cylindrical tank with a 
height-to-diameter ratio of 1.5, draft tube-to-digester diameter ratio of 0.2, draft tube-to-fluid height ratio of 0.75, 
the conical hanging baffle distance from the fluid level equal to 0.125 of the fluid height, and its outer diameter- 
to-digester diameter of 2/3. 


Keywords: Computational Fluid Dynamics (CFD), Digestion, Particles Image Velocimetry (PIV), 


Simulation 
Introduction incoming substrate and the viable bacterial 
population. These parameters are determined 
The performance of an anaerobic digester is by the flow pattern, or mixing, in the digester. 
affected by several factors, including the Complete mixing of the substrate within the 
retention time of the substrate within the digester facilitates the uniform distribution of 
digester and the degree of contact between the organisms and heat transfer. This is considered 


to be essential in high-rate anaerobic digesters 
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mixing can be performed using either 
unconfined or confined methods. In 
unconfined systems, biogas collected at the top 
of the digester is compressed and discharged 
through bottom diffusers or top-mounted 
lances (McFarland, 2001). To make the four 
gas mixing designs (Bottom diffusers, Gas lift, 
Cover mounted lances, and Bubble guns) 
comparable, MEL= 5 Wm*? at TS= 5.4% was 
used to determine the velocity of the inlet gas. 
In confined systems, the biogas is released 
through tubes. The gas lift method in a 
confined system produces the highest average 
velocity (0.080 ms‘) under the same mixing 
power (5 Wm*). In other words, mixing with 
the gas lift requires the lowest mixing power 
under the same average velocity of the flow 
field, and is the preferred method (Wu, 2014). 

The flow pattern, or mixing, inside gas- 
mixed digesters is affected by several factors 
including the biogas recycling rate, the 
clearance of the draft tube at the bottom, the 
ratio of the draft tube to tank diameter, the 
slope of the hopper bottom, the position and 
design of the biogas injection (sparger), and 
the solids loading rate (Karim ef al., 2005). 
Wei, Uijttewaal, Spanjers, Lier, & Kreuk 
(2023) assessed the impact on the treated 
sludge’s rheology as an important factor 
affecting the flow optimization and mixing 
characterization in a full-scale biogas-mixed 
digester. 

Conducting experiments to evaluate the 
effect of these parameters on mixing in the 
digester is time-consuming and _ costly. 
Therefore, simulation software like ANSYS 
Fluent is a suitable tool for designing and 
optimizing mixed gas anaerobic digesters. Wu 
(2010) presented an Eulerian multiphase flow 
model for mixing gas in digesters and 
proposed that the Shear Stress Transport (SST) 
k-w model with Low-Reynolds corrections 
would be an appropriate turbulence model to 
solve gas and non-Newtonian two-phase 
flows. 

Researchers use different indexes to assess 
the performance of their simulations and to be 
able to evaluate simulations performed with 
experimental data. Varma & Al-Dahhan 


(2007) measured the turbulence kinetic energy 
and the velocity. Karim, Thoma, & Al-Dahhan 
(2007) measured the magnitude of axial 
velocity. Wu (2010) studied the velocity 
contour, Wu (2014) used the average velocity 
and the uniformity index of velocity to 
evaluate the mixing performance, and Daplo et 
al. (2015) used the magnitude of velocity 
along the vertical axis. 

Validating the CFD simulation results is a 
necessary step. Tracer and non-invasive 
techniques are the traditional methods of 
studying gas mixing in anaerobic digesters and 
are usually used for verifying the CFD 
simulation results. Vesvikar & Al-Dahhan 
(2016), Karim et al. (2007), and Wu (2010) 
validated their models against the digester 
reported by Karim, Varma, Vesvikar, & Al- 
Dahhan (2004) and verified the flow fields 
with the measured data from Computer 
Automated Radioactive Particle Tracking 
(CARPT) and Computed Tomography (CT), a 
non-invasive technique. Dapelo, Alberini, & 
Bridgeman (2015) used Particle image 
velocimetry and a high-speed camera to 
validate an Euler-Lagrange CFD model of 
unconfined gas mixing in an_ anaerobic 
digestion. Hu ef al. (2021) proposed a novel 
approach for experimental quantification of 
mass transfer in a _ high-solid anaerobic 
digestor’s mixing process using Laser Induced 
Fluorescence (LIF) technique in a mixing tank 
equipped with multistage impellers. Flow field 
was investigated for a better illustration of the 
mass transfer, thus Particle Image Velocimetry 
(PIV) and Computational Fluid Dynamics 
(CFD) techniques were conducted for flow 
field measurement. 

The quality of mixing in a gas-lift anaerobic 
digester depends on various factors, such as 
the dimensions of the draft tube and the 
conical hanging baffle, the position of the 
baffle relative to the digester bottom, and the 
angle of the baffle. Baveli Bahmaei, 
Ajabshirchi, Abdollah poor, & Abdanan 
Mehdizadeh (2022) performed a numerical 
study and examined the influence of these 
factors on the mixing performance using 
ANSYS Fluent software. The present paper 
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extends their work by optimizing the mixing 
using the same digester configuration with 
different inlet-gas velocities. The evaluation 
criteria for optimization are average velocity, 
turbulence kinetic energy, average velocity 
gradient, and eddy viscosity of the sludge. The 
numerical results are validated using particle 
image velocimetry (PIV). 


Materials and Methods 


Methodology 

The Computational Fluid Dynamics (CFD) 
simulations were conducted using ANSYS 
Fluent software for modeling the inlet-gas 
anaerobic digester. The initial step involved 
determining the inlet-gas velocity. 
Subsequently, the effects of adding the draft 
tube and the conical hanging baffle to the 
digester design were analyzed. The 


Simulation Steps 
Determine the Inlet-gas velocity 


Simulation, Optimization and Validation Methodology 


Optimization 


Investigate the Inlet-gas Velocities 


optimization of mixing within the digester was 
achieved by varying the inlet-gas velocities 
and assessing the change in the evaluation 
indexes. The turbulence kinetic energy and the 
behavior of the sludge particles, namely their 
velocity, velocity gradient, and eddy viscosity 
were the studied indexes. The contours of the 
resulting evaluation indexes were analyzed to 
determine the optimal velocity for mixing. 
Following the simulation results, a transparent 
anaerobic digester was constructed and loaded 
with municipal sewage sludge, operating at 
optimal inlet-gas velocity. The Particle Image 
Velocimetry (PIV) method was employed to 
compare the evaluated index contours of PIV 
with those of the CFD and to validate the CFD 
simulation outcomes. A schematic 
representation of the simulation, optimization, 
and verification process is presented in Fig. 1. 


Verifying the Simualtion Results 


Determine the Conical hanging baffle's Diameter 


Determine the Height of Conical hanging baffle 


Determine the Horizontal Angle of Hanging baffle 


Determine the Length of Draft tube Investigate the Logarithmic Color contours of 


Evaluation Indexes 


Calculate the MEL 
Determine the Draft tube's Diameter Numerical Investigation of the Evaluation Indexes Construct the Simulated Digester 
i Select the Appropriate Inlet-gas Velocity 


Loading Digester by Sludge 


Perform PIV 


Fig.1. Steps used for the model simulation, optimization, and verification 


CFD simulation 


A commercial CFD software, ANSYS 
Fluent (version 19.0) was utilized to create a 
two-dimensional geometry in the design 
modeler, generate mesh, and solve the two- 
phase Eulerian model flow using the Eulerian 
multiphase approach. This two-dimensional 
model can be applied to digesters that are 
symmetrical around their vertical axis, like 
cylinders (Yang ef al., 2015). Simulations 
were performed under  unsteady-state 
conditions using Double Precision, Serial, 
Pressure-Based, and Implicit settings. The 
two-phase liquid-gas Eulerian Model of 
Viscous-SST k-omega (with sludge as the 
primary phase and biogas as the secondary 


phase) and low-Re correction were employed. 
At each time step, the iterative calculation was 
accepted as converged if all residuals fell 
below 1x10°. Final convergence was achieved 
when the average velocity of the liquid phase 
remained unchanged (Wu, 2014). 


Geometry, Computational domain, and 
mesh 

The geometry of the digester used in this 
research is based on a previously simulated 
geometry by Baveli Bahmaei et al. (2022) and 
the six steps of digester simulation are outlined 
in Fig. 1. The digester consists of a cylindrical 
tank with a flat bottom, height of 45 cm, and a 
diameter of 30 cm which results in a height to 
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diameter ratio of 1.5. The draft tube diameter 
to digester diameter is 0.2 (5 cm) and the draft 
tube height to fluid height is 0.75 (30 cm). The 
conical hanging baffle distance from the fluid 
level is equal to 0.125 of the fluid height (5 
cm), its outer diameter to digester diameter is 
2/3 (20 cm) and has a horizontal angle of 15 


Evaluation indexes 


Sludge velocity 

The velocity contour and streamlines were 
utilized in steps | to 6 of the simulation (Fig. 
1) to determine the inlet-gas velocity, draft 
tube, and conical hanging baffle 
characteristics. The uniformity of contours and 
streamlines, as well as their contribution to 
uniformity within the digester, were 
considered (please refer to Baveli Bahmaei et 
al. (2022) for more details). Sludge velocity 
was used as one of the validation indexes for 
investigating the mixing quality in a simulated 
gas-lift anaerobic digester and for selecting the 
appropriate inlet-gas velocity. The velocity 
value was compared with the sludge’s 
sedimentation velocity. Whenever the velocity 
was less than the sedimentation velocity, it 
indicated that the sludge particles would 
sediment in the digester. 


Sludge velocity gradient 
The sludge velocity gradient was used as a 


(b) 
Fig.2. The digester used for mixing optimization: (a) Geometry, and (b) Meshing; values are in cm 
and degrees 


degrees (Fig. 2). The mesh size function was 
set to curative, max face size was set to 
0.0007, and the number of nodes and elements 
were 267083 and 264281, respectively. 
Discretization error estimation was calculated 
based on the method proposed by Celik ef al. 
(2008). 


validation index for assessing the quality of 
mixing. This parameter is defined as a custom 
field function in the main menu of ANSYS 
Fluent as shown in Eq. 1 and measures the 
local velocity gradient of a mixture in 
multiphase flow using the SST k-w model as 
defined by Wu (2014). 


pwB*k 
n 


Where p and 7 are the density and the non- 
Newtonian viscosity in the liquid phase, 
respectively. 6° is 0.09 and w and k are the 
specific dissipation rate and the turbulence 
kinetic energy of the mixture, respectively. GL 
is the local velocity gradient and will be called 
the velocity gradient hereafter. 


G, 


(1) 


Turbulence kinetic energy 

Turbulence kinetic energy is used as one of 
the indexes that investigates the mixing quality 
in simulation results and is defined in Eq. 2. 
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Ou; 
Gy = —puruy = (2) 
U 


Reynolds stresses (—pu;u;) is defined in 
Eq. 3 using the Boussinesq hypothesis related 
to the mean velocity gradient. 


— Ou; Ou; 
2 
~=(pk 3 
3 (p - 


Ou; 
+ Me ix) bij 
Where 4: is the turbulent viscosity, k is the 
turbulence kinetic energy, and u (Eq. 4) is the 
velocity component. 
Uji = uU, + Uj (4) 
Where wu, and wu; are the mean and 


fluctuating velocity components respectively 
(i= 57,3); 


Sludge eddy viscosity 

Mixing quality can also be investigated 
using sludge eddy viscosity. Sludge eddy 
viscosity is the proportionality factor in 
describing the turbulent energy transfer in the 
form of moving eddies, giving rise to 
tangential stresses. Eddy viscosity is defined in 
Eq. 5 (Menter, 1993): 


_ Ur a 
vr = > @ (5) 
Mixing Energy Level 


The Mixing Energy Level (MEL) can be 
estimated using Eq. 6 (Stukenberg, Clark, 
Sandino, & Naydo, 1992). 


MEL = a 
=a (6) 


Where V denotes the effective volume of 
the digester and E denotes the energy 
consumption. Energy consumption for the gas- 
sparging (Eq. 7) was evaluated based on the 
power input formula (McFarland, 2001). 

B= P,-Q:In(=) (7) 
Py 

Where Q denotes the gas flow rate, and P 

and P2 are the absolute pressure in the tank 


headspace and at the gas-sparging inlet, 
respectively. 


Particle image velocimetry 

According to the methods used by Raffel, 
Willert, & Kompenhans (1998) and Dawkins, 
Cain, & Roberts (2012), the particle image 
velocimetry (PIV) process involves taking two 
images (I; and I) separated by time At. Both 
images were then divided into smaller regions, 
also known as sub-windows, interrogation- 
windows, or interrogation-regions. Each sub- 
window in the first image is compared with the 
corresponding sub-window in the second 
image. The sub-window with position indexes 
i and j in the first image is denoted as I,” and 
the corresponding sub-window in the second 
image is denoted as 15’. Afterward, a search 
algorithm was performed to identify a 
displacement pattern in 1;’. To do this, the 
squared Euclidean distance between the two 


sub-windows was defined in Eq. 8. 
M-1N-1 


REDS » yw (mn) 
m=0 n=0 (8) 
- iL? (m-s.n- oy 

This formula calculates the sum of the 
squared differences between all of the possible 
I,’ and I,’ sub-windows. In other words, it 
looks for the position where the sub-windows 
were the “least unlike”. Expanding the square 


parentheses in Eq. 8 would result in Eq. 9. 
M-1N-1 


R,(s.t) = > yin (s.t) 
m=0 n=0 

— iJ (m —s.n- oy 
M-1N-1 

= >, ». LJ (m. n)? (9) 
m=0 n=0 

=21 (Gb) 

-L/(m—s.n- t) 

+ LJ (m —s.n—t)? 

It should be noted that the first term, 
i (m.n)*, is a constant since it does not 
depend on s or ft. The last term, 1;/(m—s.n— 
t)*, depends on s, ¢, and only the second 
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image. To sum up, only the middle term deals 
with both of the images and this term (without 
the -2), as defined in Eq. 10, is usually referred 
to as cross-correlation (or circular cross- 


correlation). 
M-1N-1 
R(s.t) = ». > ’(m.n) 
mannee a 


-L/(m—s.n—-t) 


Results and Discussion 


The mixing conditions in the digester were 


investigated using different _—_inlet-gas 
velocities. Simulations were performed using 
inlet-gas velocities of 0.05, 0.1, 0.2, 0.3, 0.4, 
0.5, and 0.6 ms’ to study the mixing quality in 
a cylindrical digester, details of which are 
indicated in Fig. 2. 


Investigation of the evaluation indexes 

The values of the investigated indexes and 
Mixing Energy Levels (MEL) for each of the 
gas-inlet velocities are shown in Table 1. 


Table 1- Evaluated indexes and MEL for verifying the mixing quality in the gas-lift anaerobic 
digester for different inlet-gas velocities 


Turbulence Rep locit 
Inlet-gas__—_ Sludge velocity (ms"') kinetic energy highest ser y Sludge eddy viscosity (Pa s) 

velocity (m?s") eradicate) MEL 

(ms"!) min. Heed. . Gade min. a min. fies min. ag 

E-6 : E-14 . E-6 : E-17 . 
0.05 2.23 0.0236 0.30 1.0 3.8E-07 6.6 0.07 5.85 3.0E-08 0.505 
0.1 10.72 0.0291 0.43 1.0 8.3E-07 18 0.14 5.91 8.0E-07 1.01 
0.2 1.40 0.0287 0.66 1.0 64E-07 29 0.29 TAD 1.8E-05 2.02 
0.3 2.61 0.0322 0.83 8.1 0.011 359 285.23 64.50 73E-05 3.03 
0.4 3.92 0.0375 1.16 120 0.17 1398 449.68 974.12 0.65 4.04 
0.5 2.26 0.0443 1.29 4400 0.21 8370 536.97 34911.40 0.63 5.05 
0.6 9.53 0.0453 1.49 1900 0.26 5461 672.24 14858.50 0.74 6.06 
Sludge velocity Bahmaei eft al., 2022), to prevent particle 


Table 1 shows the minimum, average, and 
maximum values of sludge velocity for 
different inlet-gas velocities. The minimum 
sludge velocities were achieved in local and 
face options. The maximum velocity appears 
inside the draft tube, while the minimum value 
appears near the digester walls and at the 
bottom. The maximum velocity varies from 
0.3 to 1.49 ms’ for the studied inlet velocities 
and the average velocity only varies about 
0.022 ms"'. This indicates that the velocity of 
the particles in all internal parts of the digester 
does not increase proportionally with the 
increase in the inlet-gas velocity. This could be 
due to the formation of short-circuiting in the 
digester in areas where more mixing takes 
place. Because sludge is a non-Newtonian 
fluid and more mixing causes more decrease in 
its viscosity. 

Since the maximum sedimentation velocity 
in sludge particles is 4.7E-5 ms! (Baveli 


sedimentation, the minimum sludge velocity 
should be greater than 4.7E-5 ms’. However, 
when considering the minimum fluid 
velocities at different inlet-gas velocities, this 
goal is not achieved thoroughly at any of the 
studied inlet-gas velocities. On the other hand, 
increasing the inlet-gas velocity in gas-lift 
anaerobic digesters is limited due to the 
biological nature of anaerobic digestion. 
Therefore, a balance must be struck between 
increasing the mixing rate and reducing the 
particle sedimentation to maintain the 
conditions that prevent disruption of the 
biological process of anaerobic digestion. 


Turbulence kinetic energy 

The minimum and maximum values of 
turbulence kinetic energy for different inlet- 
gas velocities are shown in Table 1. Minimum 
turbulence kinetic energy varies between 1E- 
14 and 4.4E-11 m’s” for inlet-gas velocities of 
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0.05 and 0.5 ms‘, respectively and_ the 
maximum varies from 3.8E-7 m’s” in 0.05 
ms’! velocity to 0.26 m?s” in 0.6 ms"! velocity. 
The produced turbulence kinetic energy is very 
low for the first three inlet-gas velocities (0.05, 
0.1, and 0.2 ms‘!), has a medium value for the 
inlet-gas velocity of 0.3 ms", and is high with 


close values for the remaining three velocities 
(0.4, 0.5, and 0.6 ms”). Turbulence kinetic 
energy of different inlet-gas velocities is 
presented in Fig. 3. Higher turbulence kinetic 
energy causes more intense mixing and the 
destruction of flocs, which disrupts the 
anaerobic digestion process. 


a b Cc d e€ f g 


Fig.3. Turbulence kinetic energy contours (logarithmic color) for different inlet-gas velocities; (a) 
0.05, (b) 0.1, (c) 0.2, (d) 0.3, (e) 0.4, (f) 0.5, and (g) 0.6 ms” 


Average velocity gradient 

The average velocity gradient generates the 
turbulence kinetic energy and therefore, their 
results are similar. Results of the average 
velocity gradient for the studied inlet-gas 
velocities are presented in Fig. 4. The 
minimum average velocity gradient varies 
from 6.6E-12 to 84E-10 s"' for different inlet- 


gas velocities (Table 1). The maximum 
average velocity gradient varies from 0.07 to 
672.24 s" for inlet-gas velocities of 0.05 to 0.6 
ms". The average velocity gradient is low for 
the first three inlet-gas velocities (0.05, 0.1, 
and 0.2 ms") and high for the last three of 
them (0.4, 0.5, and 0.6 ms”). It has a medium 
value for the inlet-gas velocity of 0.3 ms". 


b Cc d e f g 


Fig.4. Average velocity gradient contours (logarithmic color) for different inlet-gas velocities; (a) 
0.05, (b) 0.1, (c) 0.2, (d) 0.3, (e) 0.4, (f) 0.5, and (g) 0.6 ms"! 


Sludge eddy viscosity 

Sludge eddy viscosity is a proportionality 
factor describing the turbulent energy transfer 
as a result of moving eddies, giving rise to 
tangential stresses. The results of sludge eddy 
viscosity for different inlet-gas velocities are 
presented in Fig. 5. The minimum and 
maximum values of sludge eddy viscosity for 
different inlet-gas velocities are shown in 
Table 1. Minimum sludge eddy viscosity 
varies from 5.85E-17 to 14.86E-14 Pa s, and 


the maximum varies from 3.0E-8 to 0.74 Pa s 
as the velocity increases from 0.05 to 0.6 ms’. 
Sludge eddy viscosity produced by the first 
four inlet-gas velocities (0.05, 0.1, 0.2, and 0.3 
ms!) has low values and the last three 
velocities (0.4, 0.5, and 0.6) are high and have 
close values. Higher eddy viscosity indicates 
higher amounts of moving eddies and high 
tangential stresses in the sludge that can lead 
to the destruction of flocs and disrupt the 
biological process of digestion. Therefore, in 
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terms of sludge eddy viscosity index, an inlet- 


gas velocity of 0.3 ms"! was appropriate. 


a b Cc d e f g 


Fig.5. Sludge eddy viscosity contours (logarithmic color) for different inlet-gas velocities; (a) 0.05, 
(b) 0.1, (c) 0.2, (d) 0.3, (e) 0.4, (f) 0.5, and (g) 0.6 ms"! 


Selecting the appropriate inlet-gas velocity 

The investigation of the evaluation indexes 
revealed that a balance between the mixing 
intensity and sludge sedimentation must be 
maintained. Higher mixing intensity can result 
in broken flocs and impairs anaerobic 
digestion. If a high inlet-gas velocity is 
selected for mixing, it can disrupt the 
biological process of anaerobic digestion. On 
the other hand, if the velocity is too low, the 
particle sedimentation rate will increase and 
proper mixing will not occur. 

Analyzing sludge velocity, turbulence 
kinetic energy, average velocity gradient, and 
eddy viscosity showed that selecting an inlet- 
gas velocity of 0.3 ms” is the most appropriate 
option. The results of CFD simulations for the 
investigated evaluation indexes for an inlet-gas 
velocity of 0.3 ms"! are shown in Fig. 6. 


a: sludge velocity 


b: turbulence kinetic energy 


The sludge velocity contour presented in 
Fig. 6 indicates that in most of the digester 
zones, zones 4 and 5 with yellow and red 
colors, the particle velocity is greater than 
1.75E-3 ms!. Considering the maximum 
sludge sedimentation velocity for the largest 
sludge particle (47 E-6 ms"! for particle size of 
2 mm) (Baveli Bahmaei ef al., 2022), particle 
sedimentation in the digester is very low. Even 
in zone 3 with a green color, the sludge 
velocity was larger than 9.9E-5 ms‘!. Only in 
zones 2 and | where sludge velocity is lower 
than 9.9E-5 ms", there is a possibility of 
sedimentation of particles larger than 0.85 
mm, which comprise 17% of the total particles 
in the sludge (Baveli Bahmaei et al., 2022). 
However, zones | and 2 cover a very small 
percentage of the digester volume, indicating 
good mixing conditions. 


average velocity gradient d: sludge eddy viscosity 


Fig.6. The resulting evaluation indexes in the digester; gas inlet velocity= 0.3 ms"! 
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Gas-sparging intensity determines — the 
amount of injected biogas for mixing and is an 
important operational assessment parameter. 
Based on the compressor’s capacity, the 
injected biogas flow rate for the inlet-gas 
velocity of 0.3 ms” was calculated to be 0.085 
m*h! in the studied digester. In the actual 
experiment, 0.085 m*h"! yielded a MEL of 3.3 
Wm”, which was close to 2.2 Wm*? that was 
reported in another full-scale gas-mixed 
digester (Dapelo & Bridgeman, 2018). 
However, this value is still much lower than 
the recommended range (5-8 Wm") needed for 
proper mixing (U. EPA, 1979). To match the 
recommended range, the inlet-gas velocity 
should be increased to over 0.7 ms”. This 
alteration requires additional investment in the 
studied digester, and the technical adjustments 
and the much higher energy consumption may 


(a) 


After selecting the inlet-gas velocity of 0.3 
ms’! as the most appropriate inlet-gas velocity, 
the particle image velocimetry (PIV) was 
performed. Due to the very dark color of the 
sludge (see Fig. 7b) and the indistinct particles 
in the images, a narrow strip of glitter was 
used along the height of the digester for PIV. 
The calculated sludge velocity, average 
velocity gradient, and sludge streamlines are 
shown in Fig. 8. The average velocity gradient 
(Fig. 8a) varies from 1.8E-6 to 34.3E-6 s', 


challenge the biogas production process 
altogether. Therefore, increasing the inlet-gas 
velocity is not an efficient strategy for 
enhancing the flow and mixing, and the 
recommended MEL criterion appears 
unsuitable for the studied digester. 


Particle image velocimetry results 

To verify the results of CFD simulations, a 
digester was constructed with transparent 
material so that photos of its inside could be 
easily taken. The transparent pilot-scale 
digester was built with the optimal 
characteristics obtained from the CFD 
simulation results and is shown in Fig. 7. It is 
made of Polymethyl methacrylate with a 
thickness of 1.5 mm. 


(b) 
Fig.7. The transparent digester: (a) empty and (b) filled with municipal wastewater sludge 


while sludge velocity (Fig .8b) varies from 0 
to 1.1x10° ms. The maximum value of 
average velocity gradient and sludge velocity 
occurred between 20 to 35 cm from the top of 
the digester, and the streamline distance is 
maximum in this zone. As shown in Fig. 8b, in 
most parts of the digester’s wall, the sludge 
velocity is greater than the minimum sludge 
velocity achieved from the — simulations, 
indicating that particle sedimentation does not 
occur. Observing the velocity contour obtained 
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from the PIV shows that the lowest velocity is 
at the junction of the wall and the bottom of 
the digester (Fig. 8b). Furthermore since there 


(a) 


are no streamlines in this area, the streamlines 


(Fig. 8c) confirm the results of CFD 
simulations. 
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Fig.8. Results of particle image velocimetry (PIV) for inlet gas velocity of 0.3 ms‘: (a) average 
velocity gradient (s"'), (b) sludge velocity (cms"!), and (c) streamline of particles in sludge 


Conclusion 


This study aimed to optimize mixing in gas- 
lift anaerobic municipal sewage sludge 
digesters. The model was built, simulated, and 
optimized, and the results were subsequently 
confirmed by building and testing the actual 
digester. 

To optimize mixing in the digester, 
different inlet-gas velocities were investigated, 
and sludge particle velocity, the gradient of 
sludge particle velocity, the turbulence kinetic 
energy, and the eddy viscosity of the sludge 
particles were evaluated. The contours of these 
evaluation indexes were analyzed to determine 
the appropriate velocity for optimal mixing, 
which was found to be 0.3 ms". 

Based on the simulation results and particle 
sedimentation velocity in the sludge, it was 
expected that the sedimentation of the particles 
would not occur in the digester at the selected 
inlet-gas velocity; except for large sludge 
particles in the small triangular section near 
the junction of the wall and the bottom of the 
digester. Subsequently, a transparent anaerobic 
digester was constructed and loaded with 
municipal sewage sludge, operating at the 
optimal inlet-gas velocity of 0.3 ms‘!. Particle 
Image Velocimetry (PIV) was employed to 


calculate sludge velocity, average sludge 
gradient, and streamlines and to validate 
simulation outcomes. According to the results 
of the Particle Image Velocimetry (PIV), in 
most parts of the digestion wall length, the 
sludge velocity is greater than the minimum 
sludge velocity achieved in the simulations. 
Moreover, the velocity contour obtained from 
the PIV shows that the lowest velocity is at the 
junction of the wall to the bottom of the 
digester and streamlines also showed that there 
are no streamlines in this area. Overall, the 
PIV method successfully validated the CFD 
simulation and showed sufficient agreement 
between the simulation and the experiment. 
The results showed that the model used for 
simulating, optimizing, and verifying the 
simulation process was successful and can be 
recommended for similar gas-lift anaerobic 
digesters, which consist of a cylindrical tank 
with a flat bottom and a height-to-diameter 
ratio of 1.5. The draft tube diameter should be 
0.2 times the digester diameter and the draft 
tube height should be 0.75 times the fluid 
height. The conical hanging baffle’s distance 
from the fluid level should be equal to 0.125 
times the fluid height, and its outer diameter 
should be 2/3 of the digester’s diameter. 
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